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We show that the process matrix in the basis of tensor products of Pauli operators or 
SU(N) generators representing low rank and sparse dynamical maps will have only a few 
distinct entries which goes as O(r^) (r is the rank). 


Understanding open quantum dynamics [1] better, developing strategies to control it [2, 3] and 
application to new technologies has gained significant interest. Quantum process tomography is 
the usual technique employed for tracking the unknown dynamics of quantum systems [4-7]. The 
prime difficulty in performing tomography tasks for large systems is the growth of resources needed 
for reconstructing the processes. For an n qubit system N = 16” — 4” independent measurements 
are required for process tomography when done in the conventional way [8]. In many situations, 
physical considerations regarding the process that is to be reconstructed and the structure of the 
problem at hand allows one to reduce the resources required to do process tomography. Indeed 
there has been considerable progress in reducing the number of independent measurements required 
for characterizing noisy processes [9-13]. In particular, extensions of the compressed sensing and 
matrix completion techniques [14, 15] that have been used effectively for quantum state tomography 
when the states are described by sparse density matrices [16, 17] to process tomography [18] provide 
substantial reduction in the resources required under certain conditions. 


I. REPRESENTATIONS OF QUANTUM PROCESSES AND PROPERTIES 

The finite time open dynamics of a quantum system is described by a quantum process (or 
operation) represented by a map, £ : p ^ £{p). Since quantum operations take density matrices 
to density matrices, the map, in general, has to be convex-linear, trace preserving and completely 
positive on all the states in its domain. The Dynamical Matrix is a matrix representation of the 
map [19] that acts on a density matrix as p^'s' £{p)rs = ^rr',ss' Pr's'- If p is d dimensional, then 
03 is a dimensional Hermitian matrix with eigenvalues Xn and eigenvectors [e^)- 

The canonical Kraus operators [20] are the matricized [21] versions of the eigenvectors of such 
that Kn = '/)fnm.sA\en)■ The mat operation stacks the elements of a column matrix row by row 
with rows of length d to generate a d x d square matrix. The number of canonical Kraus operators 
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is equal to the rank r of the dynamical matrix. The map in the Kraus (operator-sum) form is: 

r r 

^ip) = '^KnpKl where = 1 (1) 

n=l n=l 

The Kraus representation is not unique, since each can be multiplied on the left by a unitary 
matrix without violating the only constraint on the Kraus matrices given by the second equation 
in (1). Each of the Kraus operators for a map £ can be expanded in a suitable operator basis 
{Ai} as Kn = with G C. The operator basis can be chosen to be orthonormal 

(tr[74tA^] = 6jk) for convenience. Equation (1) can then be written as 

where Xij = 

ij n 

So for a given basis set {A^} the matrix x completely characterizes £. The y matrix is Hermitian 
and different Kraus representations of the same process £ have the same y matrix. From here on, 
we call the y matrix, the process matrix. 

In a typical process tomography experiment on a system made of n qubits, the {Ajj’s are 
usually taken to be the n-fold tensor products of Pauli operators. To characterize an unknown 
operation, one prepares a complete set of linearly independent input states, subject them to the 
quantum operation £ and determine the output states corresponding to each input using quantum 
state tomography. The details of performing a Standard Quantum Process Tomography (SQPT) 
can be found in [7] . We assume that the input states are initialized such that they are uncorrelated 
with the environment and hence the map that is reconstructed is completely positive. 

II. RESULT AND DISCUSSION 

We focus on dynamical maps which are sparse and low rank. Our main result is purely based 
on the properties of the process matrix. The experimental accessibility of the y matrix is the main 
motivation for studying its properties in detail for the low rank quantum operations considered 
here. Since Pauli measurements are relatively easy to perform in real quantum process tomography 
experiments, our results can be readily tested as well. For such processes we show that number of 
distinct entries in the y matrix goes as O(r^). 

The y matrix can be written as an outer product as follows 

r 

X = Y.£A (3) 

i=l 

where 

= [Kf \ Kf \ ..., with Kf) = tr(KiArf), 

and 

2i = [Y\Y\ ..., Y ], with Y = ti{YXd). 
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Here {Aj} denotes the convenient operator basis used for defining the x niatrix and the superscript 
T denotes the transpose operation. For the discussion that follows we take the basis to be made 
of the n-fold tensor products of Pauli matrices. 

The product is Hermitian and so is x- The elements of are conjugates of Tj. Let us 
assume that the Kraus matrices, K,, each have only a maximum of r non-zero entries. For an 
s-sparse matrix, it’s rank satisfies the bound 0(1) < r < s and also a tighter bound 0(1) < r < 
min(s,(i) for a d dimensional matrix. [22]. The Tj, which is the the matrix constructed out of the 
trace with X/s of the Kraus matrices, will have 0(r) distinct entries. Consider the trace of an 
arbitrary 2x2 matrix with the Pauli matrix ui 


/ a (5)\ /O l\ 

\® y ) 0) 


b + X 


Only two elements b and x (the circled ones) contribute to the trace. Since the trace operation 
is with Pauli matrices or their tensor products, it is evident that only d elements in the case of a 
dxd matrix will contribute to the trace, since the Pauli matrices or their tensor products have only 
d non zero entries. Consider an r-sparse Kraus matrix which means that it has only r non-zero 
entries. 


The probability of an entry being non- zero = 

d^ 

The probability that only one element contributes to trace 

r \'^-i r 


= ^Ci 


1 


d^J 


The probability that only one element from each of the r Kraus operators contributing to the 
trace is high. Remember that r « d^ which validates the high probability argument. If they are 
located in such a way that only one element contributes towards its trace with a particular Aj, then 
4r different possibilities alone arise, since the trace with Pauli tensors produces ±1 or ±i alone. 
This indicates that if more elements contribute to the trace with a particular Aj, the number of 
distinct entries scales linearly with r. 

From Eq. (3), it is therefore clear that the number of distinct entries in the x matrix is additive 
w.r.t to the Kraus matrices and since they are r in number, the number goes as O(r^). It can be 
understood from this construction that if different Kraus matrices have the same matrix positions 
of non-zero entries the number of distinct entries will not increase. 

Consider the case where only a single non zero element is present in the Kraus operators, but 
their positions are different. 


Ml 



M2 = 



^ af 0 0 

0 02 —0 
0 ia2 02 0 

^ o^ 0 0 of j 


XI = 






The X matrix has two non zero entries. 
Consider a second case as follows 
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Li = 




0 

02 



+ “i 

0 

0 

a\ — 02^ 


0 

0 

0 

0 


0 

0 

0 

0 


- ai 

0 

0 

al + al ) 


The X matrix has two non zero entries. 

It is worth noting that the x matrix in the basis of SU(N) generators will also have the same 
nnmber of distinct elements due to the fact that the trace of Kraus matrices with generators of 
SU(N) Lie algebra also produces ±1 or alone.This means that the x matrix written in the basis 
of SU(N) generators will also have only O(r^) distinct entries. 


III. NUMERICAL SIMULATION 

The histograms shown below clearly indicates that the number of distinct entries in the x matrix 
are peaked at low numbers which is less than r^. The histograms were made by constructing x 
matrices from sparse Kraus matrices, stemming from sparse and low rank dynamical matrices 
and counting the number of distinct entries in absolute values. For rank 3, we created three 8x8 
matrices each having only 3 entries by allocating the non zero entries in any three random positions 
of each of the three matrices such that Eq. 1 is satisfied and the number of distinct elements in 
absolute values was counted for the x matrix developed in the basis of CTj (8) Oj ® (Tk- This was 
repeated for 100,000 realizations by choosing all possible permutations of the location of the non 
zero entries which satisfied Eq. 1. A similar process was repeated for rank 2 case as well where we 
have only two Kraus matrices to start off with. 



FIG. 1. (Color online) Histogram of 100,000 realizations showing the number of distinct matrix elements 
(in absolute values) of a 64 dimensional x matrix corresponding to a dynamical matrix of rank 2. 
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FIG. 2. (Color online) Histogram of 100,000 realizations showing the number of distinct matrix elements 
(in absolute values) of a 64 dimensional x matrix corresponding to a dynamical matrix of rank 3. 

IV. CONCLUSION 

We have shown that when the underlying process is represented by a sparse dynamical matrix, 
then in the x matrix form, a lot of entries are repeated. Thus a x matrix reconstructed from process 
tomography data with many repeated entries is another indicator for the underlying process to be 
low rank and sparse. 
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